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We use the unitary and analytic model of the pion electromagnetic form factor in order to evaluate in the lowest 
order the e + e~ — > tt + tt~ contribution into the muon magnetic anomaly. We demonstrate, that this technique 
enables us to reduce the uncertainty of the theoretical prediction importantly in comparison to usual approaches, 
where the measured data are integrated directly. 



1. Introduction and Motivation 

The question of an experimental evidence of 
physics beyond the Standard Model (SM) has 
been and still is an issue in particle physics. This 
topic is especially popular nowadays, when new 
high-energy experiments are being constructed 
and put in work. Regarding the SM as a low- 
energy approximation of a more general theory, 
one can expect the new physics to manifest itself 
in two ways: important effects in measurements 
at high energies or small effects at relatively low 
energies. The latter is the case of the muon mag- 
netic anomaly, often referred to as the muon g — 2 
problem. To answer the question, whether some 
new physic plays a role in this case, one needs 
to achieve a high precision in both, theoretical 
prediction and experimental measurement. 

The muon gyromagnetic ratio relates the 
muon spin and magnetic momentum ~jt = 

(e/2m M ) if and the magnetic anomaly is de- 
fined by a M = — 2) /2. The higher muon 
mass makes its magnetic momentum, in compar- 
ison to the magnetic momentum of the electron, 
more sensitive to the non-perturbative hadronic 
effects and thus can the determination serve 
as a good test field of the theory of strong in- 
teractions. A high-precision experimental result 
is available, the anomaly et M was measured by 
the g — 2 collaboration in the E821 experiment 
at BNL with unprecedent precision of 0.7 ppm 
PQ. A high-precision theoretical prediction is 



available also, but only in the electromagnetic 
and electroweak sector the accuracy is satisfac- 
tory. In case of the hadronic contribution a'^ ad 
more precision is desirable, because this contri- 
bution is the dominant source of the total uncer- 
tainty of the theoretical prediction, which does 
not allow to conclude on possible physics beyond 
the SM. Depending on the author and the an- 
alyzed data one can cite different numbers. If 
one takes the values from the summary publi- 
cation [2] af p = 116 592 093(63) x 10" 11 and 
af = 116 591810(210) x 10~ n one observes 
Aa„ « 1.3a A(V , Aa M = \af - a»*|. 

The hadronic part is dominated by the lowest 
order and this can be calculated using the formula 

<a , LO iray r *jl k{s)r{s) , (d 

3 V 7T/ J 4m 2 s 

where K(s) = £ dx[x 2 (l - x)]/[x 2 + (1 - x)^] 

and R(s) = a LO (e+e" -> had) /^f. The 
most important contribution comes from the pion 
channel a^ l ad ' LO (e + e~ — » tt + tt~), which we study 
in this work. The data are used to evaluate the 
integral up to few GeV, above a perturbative cal- 
culation is possible. In order to compare to other 
authors we chose three different limits for the up- 
per integration limit (3.24GeV 2 , 2.0449GeV 2 and 
0.8 GeV 2 ). 

The expression ([T]) was up to now [314] always 
evaluated only by the direct integration of the 
clo (7r + 7r - ) data (with small exceptions, see [S]). 
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Our motivation and our aim is to demonstrate, 
that the use of an appropriate model allows for 
a dramatic error reduction. The error reduction 
in this domain is of a crucial importance, since 
there is an effort to decrease even more the ex- 
perimental error [6j and so, if one wants to solve 
the g— 2 puzzle, the theoretical precision will have 
to follow. 

2. Unitary and Analytic Model of the Pion 
Electromagnetic Form Factor 



We base our evaluation of 



the unitary and analytic (U&A) model of the 
pion electromagnetic (EM) form factor and relate 
the form factor to the cross section er tot (e + e~ — > 

7T+7T-) = 2i£ (1-Aml/ty \F n (t)\ 2 , so as to ob- 
tain R(s). For a reliable prediction one needs 
to construct a model that incorporates all known 
properties of the pion EM form factor. The form 
factor is charge-normalized F w (0) = 1 and pQCD 
predicts its asymptotic behavior -FV(£)|t|-+oo ~ 
t -1 , where t = q 2 = —Q 2 - It is known, that 
F„(t) is an analytic function in the whole com- 
plex t plane, besides the branch points and the 
cut on the positive real axis, the cut going from 
the lowest branch point t = 4m^ up to +oo. 
The pion form factor has also cut on the sec- 
ond, unphysical Riemann sheet for — oo < t < 
and satisfies the unitarity condition Im [F„(t)] = 
[A{(t)]* F n (t) + a(t), where A\(t) is the P -wave 
isovector transition amplitude for elastic ir + n~ 
scattering and a(t) approximates all higher con- 
tributions. For 4m 2 < t < {m v + m u ) 2 one has 
a(t) =0 and arrives to the so-called elastic uni- 
tarity condition. 

The construction of the U&A model starts from 
the Vector Meson Dominance (VMD) picture in 
which 
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fv 



v=P,p',p" 

where f V7r7r is the vector meson-pion coupling con- 
stant, /„ is the universal vector meson coupling 
constant and three resonances are taken into ac- 
count. We proceed to a non- linear transformation 

4(t in - t a ) 



t = ta- 



in order to build in the cut, the lowest branch 
point to and an effective branch point £j n , the 
latter meant to approximate all higher branch 
points. The model then becomes 
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v—P,p',p" 



(W N - W vQ ) {W N + W v0 ) 
(W - W v0 ) (W + w v0 ) 



{W N - 1/W v0 ) (W N + 1/Wvo) 
(W - 1/Wvo) (W + 1/W v0 ) 



{fvTnrffv] 



where W N = W(t)\ t=Q and W v0 = W(t)\ t=m 2. 
The expression can be modified further; it can be 
shown that depending on the relative positions of 
thresholds and masses one has to < m v < hn => 
W v = -W* and U n <m v ^ W v = l/W*. The 
masses and U n are considered as free parameters 
of the model and so their relative positions are 
not fixed. However, in order to provide explicit 
formulas in the following text, we will here sup- 
pose to < m p < ti n < mp 1 , Trip" , what is actually 
being confirmed by the result of the data fitting. 

The next step in the model construction is in- 
corporation of the cut on the second Riemann 
sheet. We use a succession of poles and ze- 
ros generated by a rational function to approx- 
imate this cut (Pade approximation). We do it 
by adding to the model a multiplicative term 
^-wi^w-Zr) > where Wz and W P are free pa- 
rameters from the interval < Wz,p < 1 (corre- 
sponds to — oo < t < 0). The construction of 
the model is finally achieved by giving the vec- 
tor mesons non-zero decay widths W v o — > W v = 
W{t)\ / -iv \2. The model gets form 

F mm ( i-w 2 \ 2 {w-w z )(w N -w P ) 
W(*)J - ( t-fjtH {Wn _ Wz){w _ Wp) 



i 



w 2 N 



{w N - w p ) (w N - w;) 
(W - W p ) (w - w*) 

{w N - i/Wp) (w N - i/w;) 



[l/W-W] 2 



(W- 

E 



- l/Wp) (w - l/W*, 

(W N - W v ) (Wn - W*) 
(W-W V )(W-W*) 



{fpTnr/ fp) 
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3. Error Evaluation and Results 
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Figure 1. Fit of the data with the U&A model. 
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The ratio of couplings fpmr/fp is a free parameter 
of the model, the two remaining ratios, fp'w/fp 
and /p"7T7r//p", can be related to the first one 
by using the normalization condition i^(0) = 1 
and by considering the behavior of the imaginary 
part of the form factor at q = 0. The p — uj 
interference is taken into account when fitting the 
experimental data by a Breit-Wigner term, the 
fitting function is 

F„ [W(t)} + Ref- 



ill- 



t 



where d> = arctan • 



■i p p i and the amplitude R is 

,.ip— m£ 

considered as an additional free parameter. The 
fitting function thus has 11 free parameters in to- 
tal (U n , three masses, three widths, ratio fpmr/fp, 
W z , Wp and R) and is used to fit 523 experi- 
mental points measured in experiments CLEO, 
NA7, OLYA, CMD, CMD-2, SND, KLOE, in 
JINR Dubna and in F T collaboration at Jefferson 
Lab 17l8l9ll0lllll2ll3ll4ll5ll6ll7ll8ll9l20l2ll^2l 
I23l24l25l26l27l28l29l30l3ll32l33l34| . The best fit 
result is shown in Figure [TJ 



Two approaches were used for the error eval- 
uation. The first one was a "Monte Carlo" 
method based on a random number generator 
with the assumption of the Gaussian distribu- 
tion for the uncertainties of the published data 
points. For each point a new one was ran- 
domly generated using the Gaussian probability 
density function with the mean identical to the 
original point and a equal to the published er- 
ror. Doing this for each point, new "random" 
data set was obtained. This data set was then 
fitted and the values of the parameters of the 
model pi as well as the value of a^f' d,LO (tt + tt~) 
were extracted. Repeating the whole procedure 
4000 times, we reached statistics high enough to 
allow us for a reliable error calculation. The 
mean ap iad ' LO (tt + it~) and the er were calculated 
from the 4000 values and since the mean is 
not, in general, identical with the optimal-fit 
value a h a a QPT ( 7r+7r ) we present asymmetric un- 



certainties a 



ted , i O ( ^- ) = a WO (7r+7r _ ) -A 



where A = a 



had.LO / 4- — \ 



had,LO 



(7T + 7T ) 



had,LO / 4- — \ 

■ a u (tt^tt ) - a 



had,LO i 4- — \ 
yTT 7T ). 



and B — w -p u, n yu n f u,^ QpT 

In the second approach the program MINUIT 
was used to establish the uncertainties of the 
model parameters. Then, taking the numeri- 
cal derivatives for ■^-a p iad LO (tt + 7t~), the uncer- 
tainty was propagated to a ^ ad ' LO (n + Tr~) using 
the covariance matrix. In this method the errors 
are symmetric. 

In addition to the integration of the model, we 
also performed a direct integration of the data 
points based on the trapezoidal rule, so as to 
cross-check our compatibility with other authors. 
Our results and some results from other authors 
|3|4|5j are summarized in Table Q] 

4. Discussion, Summary and Outlook 

The use of the U&A model dramatically re- 



duces the error on a l }, ad ' LO i 



). This is not 



an arbitrary feature of the model but originates 
from model-independent information which is ad- 
ditional to the data in the integration region and 
which can be taken into account when the model 
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a j*.d,LU( e + e - — „.+„--) x 10 n 



Interval \GeV 2 ] Ami < t < 3-24 Ami < t < 2.0449 Ami <t<0.8 

Model, method 1 5132.36+g 5128.22+g;^ 4870.24+g;|g 

Model, method 2 5132.37 ± 3.00 5128.25 ± 2.86 4870.44 ±2.64 

Data 5035.33+^32 5031 - 22 +28:94 ^ G - 77 +27.li 

Davier 5040.00 ± 31.05 

Hagiwara et al. 5008.2 ± 28.70 

Yndurdin et al. 4715 ± 33.53 

Table 1 

Our results and results of other authors [31415] . 



is used. The most important sources contributing 
to error reduction are 

• Expected smoothness of the F n (t) at small 
scale At: The model provides a function 
behaving smoothly at small At. 

• Experimental data outside the integration 
region: The fit is done not only to the data 
inside, but also to the data outside the in- 
tegration region. 

• Theoretical knowledge on F w (t): The model 
respects all known properties of the pion 
EM form factor. 

Especially the first point plays an important role. 
The new precise data tend to lie above older, less 
precise data and, in some regions, the vertical 
spread of the data is very important, at the limit 
of inconsistency. If the calculation of the inte- 
gral is based directly on data, then less precise 
data shift the mean value of the integral and en- 
large the uncertainty. When the model is used, 
the predicted behavior of F„(t) as given by the 
result of the fit is mostly determined by precisely 
measured points and is only little influenced by 
data with important uncertainties. This leads to 
more appropriate mean value and smaller errors. 
In consequence, one arrives to the mean value of 
a had,LO ^ 7r + n --j which j s higher then what is ob- 
tained by the direct data integration and to much 
reduced uncertainty. The shift in the mean value 
goes in the "right" direction and brings the the- 
oretical value closer to the experimental one. 



The error estimates from the two used methods 
are not fully compatible, the first method gives 
smaller errors. This might be related to statis- 
tical fluctuations (1st method) and to approxi- 
mations - numerical derivatives and linearization 
(2nd method). 

In this article we presented the calculation of 
a had,LO (tt+ti--) b asec i on the U&A model. This 
approach allows for important error reduction 
and we plan to use it also for evaluating the con- 
tributions of other channels to a t }, ad • . 
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